BIOMETRIC STUDY OF 
BARREMITES SUBDIFFICILIS (KARAKASCH) 

by R. A. REYMENT and P. SANDBERG 


Abstract. A study is made of allometric growth in four dimensions of the ammonite Barremites subdifficilis 
(Karakasch) from the Neocomian of the Crimea. The four dimensions are further analysed simultaneously by 
the multivariate statistical method of principal component analysis. The results are compared with those found 
for two species of Triassic ammonites. 

It is often assumed by quantitative paleontologists that there will be an allometric 
growth relationship between many of the dimensions of the coiled ammonite shell, 
Obata (1959) investigated several species of Desmoceratidae and found almost isometric 
growth for most combinations of the variables studied but in a certain number of cases 
possibly significant positive and negative allometric growth. 

In the present analysis the variables studied are maximum diameter of shell (%), 
maximum diameter of umbilicus (x 2 ), height of last whorl („y 3 ), breadth of last whorl („y 4 ); 
all measurements are in millimetres. The measurements were converted to base 10 
logarithm values for the bivariate and quadrivariate studies. 

The material was collected by Dr. V. V. Drushshits, Department of Paleontology, 
University of Moscow, who also checked the measurements. The specimens are stored 
in the University of Moscow; all (32) specimens were derived from the same bed in the 
Barremian sequence of the Crimea, U.S.S.R. The quadrivariate computations were 
performed on the computers of the Swedish Board for Computing Machinery, who 
provided free machine time. Part of the study was supported by a grant from the Swedish 
Natural Science Research Council. This study was completed during the junior author’s 
tenure of a fellowship from the Southern Fellowships Fund, Chapel Hill, North 
Carolina, U.S.A. 


RELATIVE GROWTH ANALYSIS 

Standard least squares regression was used in the allometric analyses. The equations 
between pairs of variables are: 


x x = 7-04x 2 076 

(1) 

x x = 2-43xr 3 0 ' 93 

(2) 

Ay = 2-99 x 4 °" 

(3) 

x 2 = 0-39 a- 3 104 

(4) 

x 2 = 0-45V 14 

(5) 

x 3 = 1-43 x 4 °" 

(6) 


The exponents of equations (2), (3), (4), and (6) differ so slightly from unity that there 
can be little doubt that growth in these cases is isometric. Equations (1) and (5) are 
indicative of possible allometric growth. The 95 per cent, confidence intervals for the 
regression coefficients of these two equations were computed. For equation (1) it was 
found that 0-69 < /? < 0-84, where /3 is the population regression coefficient. Thus there 
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is only a 5 per cent, chance of the population coefficient lying outside of these limits; 
hence it may be concluded that growth here is of the negative allometric kind. For 
equation (5) the 5 per cent, confidence intervals for the regression coefficient are 
1*07 < j8 < 2-00. The lower limit does not include unity and it may be suggested that 
a slight positive allometric relationship exists between x 2 and y 4 . 


PRINCIPAL COMPONENT ANALYSIS 


The covariance matrix for the logarithmically transformed observations of the four 
variables is (N = 32) 


S = (32-1)- 1 


0-5942 

0-6753 

0-5975 

0-5547 

0-6753 

0-8831 

0-6631 

0-6393 

0-8831 

0-6631 

0-6405 

0-5573 

0-5547 

0-6393 

0-5573 

0-5626 


The eigenvalues and eigenvectors of this matrix were extracted with the aid of an elec¬ 
tronic computer. They are given in the following table. 



Percentage of total 

Eigenvalues 

variation 

A 4 = 0-08166 

94-44 

A 2 = 000314 

3-63 

A 3 = 000131 

1-50 

A 4 = 0-00037 

0-37 


Eigenvectors 

b\ - (0-480, 0-570, 0-486, 0-458) 
b' 2 - (0-194, -0-790, 0-547, 0-200) 
b' 3 = (-0 065, -0-169, -0-522, 0-833) 
b' 4 = (0-853, -0-153, -0-438, -0-239) 


The first principal component accounts for almost all of the variation. Together with 
the second component, 98 per cent, of the total variation may be explained. It is interest¬ 
ing to observe that all variables are roughly equally balanced in the first principal com¬ 
ponent, which implies that the four dimensions participate equally in the variability 
of the shell morphology. The second principal component is mainly concerned with 
variation in log x 2 and log x 3 ; that is there is a shape relationship between umbilical 
diameter and height of whorl. The third principal component is representative of varia¬ 
tion in height and width of whorl; it could be shown to be statistically significant by 
Lawley’s test (1956). 

To the extent of the writers’ knowledge, this is the first time a principal component 
analysis has been made on the test of an ammonite and there is therefore no compara¬ 
tive literature available. In order to gain some impression of the importance of the above 
figures, two principal component analyses were run on Triassic ammonites, the data 
being taken from Kummel and Steele (1962). The observations were transformed 
logarithmically as before. 

(a) Dieneroceras spathi Kummel and Steele (Kummel and Steele 1962, p. 659). The 
material (N = 42) derives from Elko County, Nevada, and Cottonwood Canyon, Cali¬ 
fornia (zone of Meekoceras grcicilitatum). 

The covariance matrix is: 


S = (42—l)- 1 


(0-5419 

0-5831 

0-3770 

0-4764 

0-5831 

0-7772 

0-4007 

0-4747 

0-3770 

0-4007 

0-3058 

0-3490 

,0-4764 

0-4747 

0-3490 

0-4738 



R. A. REYMENT AND P. SANDBERG: BIOMETRIC STUDY OF BARREMITES 729 


The eigenvalues and eigenvectors of this matrix are presented in the following table. 



Percentage of total 

Eigenvalues 

variation 

Aj. = 0-046445 

90-73 

A 2 = 0003614 

7-06 

A 3 = 0-000729 

1-42 

A 4 = 0-000399 

0-78 


Eigenvectors 

b\ = (0-526, 0-603, 0-377, 0-467) 
b 2 = (0-153, -0-754, 0-282, 0-573) 
b' 3 - (-0-241, -0-027, 0-881, -0-405) 
b 4 = (0-801, -0-260, -0-036, -0-538) 


(b) Meekoceras gracilitatum White (Kummel and Steele 1962, p. 694); the material 
is from Elko County, Nevada, U.S.A. (Zone of M. gracilitatum ); it comes from several 
beds. The covariance matrix is (N — 43): 


S = (43-1)- 1 


/T-1672 
1-5413 
0-9800 
[ 1-0675 


1-5413 

0-9800 

1-0675 

2-2879 

1-3173 

1-3681 

1-3173 

0-8604 

0-9112 

1-3681 

0-9112 

1-0649 


The eigenvalues and eigenvectors of this matrix are given below. 


Percentage of total 
Eigenvalues variation 


Eigenvectors 


A x = 0-121998 

95-23 

A 2 = 0-004907 

3-74 

A 3 = 0-000799 

0-62 

A 4 - 0-000517 

0-40 


b\ = (0-471, 0*653, 0-403, 0-434) 
b' 2 = (0-221, -0-706, 0-187, 0-647) 
b' 3 = (0-475, -0-273, 0-559, -0-623) 
b 4 = (-0-709, 0-031, 0-700, 0-074) 


It will be observed that in both cases the first principal component accounts for more 
than 90 per cent, of the total variation. Moreover, the contributions of each of the dimen¬ 
sions are approximately equal; this was also found to be so for Barremites subdifficilis. 
In fact, if one computes the angle between the theoretical vector k' = (0*5, 0-5, 0-5, 0-5) 
and each of the first eigenvectors for the three species here discussed, it is found that it is 
in all cases roughly zero degrees, which implies accordance with the hypothesis of equal 
contribution by all dimensions. 

The second principal component is for the Cretaceous ammonite mainly indicative 
of covariation in variables log x 2 and log x 3 , while for the two Triassic ammonites 
log x 2 and log x 4 covary. For all three species the third principal component represents 
covariation in variables log x 3 and log x>. 

These results are enlightening. They disclose the perhaps surprising situation that 
diameter of shell, umbilical diameter, height of whorl, and breadth of whorl are possibly 
equally important in the morphologic variation of ammonites. The remaining principal 
components tend to differ with respect to the sign and magnitude of the components of 
the eigenvectors connected with them. It is conceivable that these components may be 
of specific significance and that the vector of principal components could be of taxo¬ 
nomic use in characterizing a species. In the following table we supply the mean vector 
for the principal components of each species (computed from the equation y = B'z, 
where B is the matrix of eigenvectors and (z 1? .... z 4 ) = (log x l5 ..., log x 4 )). 


Species 

Barremites subdifficilis 
Meekoceras gracilitatum 
Dieneroceras spathi 


Mean vector of principal components 

y x = (1-974, -0-316, 0-075, 0-351) 
y 2 = (2-433, 0-920, 0-314, -0-310) 
y 3 = (2-127, 0-317, -0-056, 0-323) 
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CONCLUSIONS 

The growth relationships between the variables maximum diameter of shell, umbilical 
width, height of last whorl, and breadth of last whorl are for Barremites subdifficilis 
(Karakasch) almost isometric. Slight negative allometric growth could be shown to 
exist between maximum diameter and umbilical diameter and slight positive allometric 
growth was shown to occur between umbilical diameter and breadth of the last whorl. 

The covariance matrix of these variables was subjected to the multivariate statistical 
method of principal component analysis, whereby the total variation is broken up into 
components of variation. In order to provide our conclusions with a somewhat more 
general basis two species of Triassic ammonites were also analysed by the same method. 
The results show that almost all of the biologic variation in the shells of these ammonites 
with respect to the variables studied is contained in the covariation of these variables 
studied and that each of these contributes equally to the variation. This result would 
appear to vindicate the common use of these four dimensions in the description of 
ammonites. It would seem possible that the mean vector of principal components might 
be useful in the characterization of species. 
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